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OF  THE  ORBIT  OF  COSMOS  837  ROCKET 


I  C)  H. I  Hiller 


SUMMARY 


The  orbit  of  Cosmos  837  rocket  (1976-62E)  has  been  determined  at  36  epochs 
between  January  and  September  1978,  using  the  RAE  orbit  refinement  program  PROP6 
with  about  3000  observations.  The  inclination  was  62,7°  and  the  eccentricity 
0.039,  The  orbital  accuracy  achieved  was  between  30  m  and  150  m,  both  radial  and 
cross  track. 

The  orbit  was  near  29:2  resonance  in  1978  (exact  resonance  occurred  on 
May  14)  and  the  values  of  orbital  inclination  obtained  have  been  analysed  to 
derive  lumped  29th-order  geopotential  harmonic  coefficients,  namely: 


g-0.2 

10  C29 


10  ±  15 


o_0,2 

i°  s29 


-  -  76  ±  12  . 


These  will  be  used  in  future,  when  enough  results  at  different  inclinations  have 
accumulated,  to  determine  individual  coefficients  of  order  29,  The  values  of 
lumped  harmonics  obtained  from  analysis  of  the  values  of  eccentricity  were  not 
well  defined,  because  of  the  high  correlations  between  them  and  the  errors  in 
removing  the  very  large  perturbation  (31  kra)  due  to  odd  zonal  harmonics. 
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1  INTRODUCTION 

Cosmos  837  rocket,  1976-62E,  entered  the  following  orbit  on  1976  July  I:  inclination 
62.75°,  perigee  height  440  km,  apogee  height  920  km,  period  98.4  min  and  eccentricity 
0.034.  Decay  is  expected  about  1984. 

The  orbit  seemed  promising  for  analysis  to  obtain  29th-order  'lumped1  geopotential 
harmonic  coefficients,  by  virtue  of  the  relatively  slow  passage  through  29:2  resonance 
with  the  Earth's  gravitational  field  -  when  the  satellite's  track  over  the  Earth's  surface 
repeats  every  two  days  while  the  satellite  makes  29  revolutions  of  the  Earth.  The  orbit 
has  been  analysed  for  the  nine-month  period  1978  January  to  September  when  the  effects  of 
29:2  resonance  were  significant:  exact  29:2  resonance  was  on  1978  May  14. 

The  orbit  was  determined  at  36  epochs  using  the  RAE  computer  program  PR0P6 

The  inclination  and  eccentricity  values,  after  clearance  of  perturbations,  were  fitted 

7 

by  least-squares  theoretical  curves,  using  the  THROE  program^.  The  best  fittings  gave 

lumped  29th-order  coefficients  which  are  nominally  the  most  accurate  so  far  obtained. 

These  lumped  coefficients,  with  others  at  different  inclinations,  will  be  used  to  obtain 

individual  29th-order  harmonic  coefficients.  A  few  other  values  have  previously  been 
3~ 5 

determined  ,  but  several  more  are  required,  especially  at  inclinations  near  90°. 

2  ORBIT  DETERMINATION 

2. 1  Observations 

Over  3100  observations  were  available  for  the  36  orbits  selected  for  determination, 
over  the  period  1978  January-Sep tember .  \bout  15%  of  these  observations  were  rejected 
due  to  not  fitting  well,  to  leave  a  working  average  of  74  observations  per  orbit.  Three 
of  the  orbits  had  the  benefit  of  highly-accura te  Hewitt  camera  observations. 

The  largest  group  of  observations  used  was  about  2300  from  the  US  Navy;  a  further 
400  came  from  British  radar,  130  from  the  kine theodoli te  at  the  South  African  Astronomical 
Observatory  (SAAO) ,  in  the  southern  hemisphere,  32  from  the  theodolite  in  Jokioinen, 
Finland  and  nearly  300  observations  were  from  volunteer  visual  observers,  supplied  by  the 
Anpleton  Laboratory  at  Slough. 

In  the  course  of  the  orbit  determinations,  it  was  noticed  that  several  groups  of  the 
South  African  observations  were  rejected,  many  with  consistent  errors  of  -3  or  -4  seconds. 
Correction  of  these  timing  errors  led  to  satisfactory  acceptance;  subsequently,  the  SAAO 
confirmed  that  such  errors  had  been  made  in  these  observations  (and  possibly  others)  over 
a  5-week  period  at  the  beginning  of  1978. 

2 . 2  The  orbits  and  their  accuracy 

The  36  computed  orbits  are  given  in  Table  1,  where  it  can  be  seen  that  the  sd  in 
inclination,  i  ,  varies  from  0.0003°  to  0.0017°,  the  rms  value  being  0.0008°;  by  compari¬ 
son,  the  rms  for  the  three  Hewitt  camera  runs  is  0.0005°.  For  eccentricity,  e  ,  the  sd 
— 6  —6 

varies  from  4  *  10  to  19  y  10  .  The  best  sd  values  of  i  and  e  are  equivalent  to 

30  m  in  position.  For  the  right  ascension  of  the  node,  II  ,  the  average  sd  is  0.0009°, 
equivalent  to  100  m. 


Table  1 

Orbital  parameters  for  Cosmos  837  rocket,  with  standard  deviations 


r  1 - 

l  MJD 

Date 

1978 

a 

e 

i 

- 

M0 

M2 

0 

l 

1  S 

l  14)514.0 

Jan  6 

7056.6206 

0.038256 

62.7406 

15.7015 

2.181 

14.924 

52  7  3.0814 

0.02)74 

0.52 

8.9 

61 

2 

14 

1 1 

9 

10 

10 

3 

10 

2  522.0 

14 

7056.3017 

0,038355 

62.7417 

349.9948 

'3.510 

81.071 

527  3.4390 

0.01961 

0.44 

H.9 

6  3 

2 

8 

8 

8  j 

8 

8 

2 

8  | 

3  5)2.0 

24 

7055.9888 

0.038501 

i  62.7420 

317.8549  1 

5.184 

1  257.305 

5273.7898 

0.01666  1 

0.56 

8.6 

76 

3 

7 

B 

10 

10 

10 

3 

10  ( 

4*  540.0 

Feb  1 

7055.6339 

0.038559 

62.7423 

292.1450 

6 . 504 

329.042 

5274.1879 

0.0318)  I 

0.50 

7  ,h 

65 

3 

8 

3 

8 

j  10 

10 

3  i 

16 

5  548.0 

9 

7055.2032 

0.038676 

62.7386 

266.4234 

1  7.818 

44.588 

5274.6707 

0.02684 

0.49 

7.6 

64 

555.0  16  7054.9588  0.038812  62.7382 

3  6  I  6 

564.0  25  7054.7220  0.038966  62.7377 

3  19  17 

574.0  Mar  7  7054.3304  0.039081  62.7388 

3  13  12 

584.0  17  7053.8633  0.039178  62.7431 

2  5  8 

594,0  27  7053.4390  0.039300  62.7375 

3  8  13 

604.0  Apr  6  7052.7617  0.039338  62.7345 

4  9  8 

611.0  13  7052.0397  0.039385  62.7365 

*  5  5 

616.0  18  705 1 . j292  0.039357  62.7385 

6  119 

624.0  26  7050.8419  0.039424  62.7402 

2  119 

633. 0  May  5  7050.1571  0.039559  62.7371 

!  6  17  U 

j  641.0  13  7049.6910  0.039619  62.7373 

2  7  6 

j  648.0  20  7049.3633  0.039675  62.73S6 

f  10  7 

!  655.0  27  7049.0539  0.039729  62.7389 

|  1  4  S 

i  665.0  Jun  6  7048.5968  0.039820  62.7417 

I  3  13  11 

J  673.0  14  7048.3420  0.U39880  62,7421 

j  18  6 

681.0  22  7047.9572  0.039976  62.7381 

I  5  12  8 

I  688.0  29  7047.3403  0.040027  62.7362 

1  9  5 

695.0  Jul  6  7046,7772  0.040034  62.7372 

2  15  9 

701.0  12  7046.3310  0.040034  62.7387 

3  14  6 

707.0  18  7045. 770J  0.040048  62.7407 

4  16  10 

!  714.0  25  7045.1294  0,040104  62.7388 

4  1)  8 

I  710.0  29  7044,8866  0.0401)8  62.7  384 

2  7  5 
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47.5931 

9 
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133.3921 
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1  14.0078 

16 

91 . 385  ) 
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78.4571 


l  8.947  248.384  5274.9447  0.01785  0.61  6.9 

7  j  7  3  14  | 

10.464  204.080  15275.2103  0.01677  0.64  9.3 

17  16  j  4  15 

12.042  38.375  5275.6499  0.02384  0.53  8.4 

11  11  3  14 

13.683  237.579  5276.1743  0.02422  0.60  9.1 

8  8  3  6 

15.405  81.595  5276.6500  0.02807  0.77  9.4 

15  14  3  13 

16.990  291.491  5277.4101  0.05014  0.54  5.9 

10  9  4  21 

18.072  156.1)0  5278.2209  I  0.06162  0.6?  5.1 

8  8  4  20  j 

18.934  268.728  5278.7942  0.05279  i  0.49  5.7 

12  tl  7  (  28 

20.228  22.365  5279.5664  0.04268  0.47  8.9 

11  10  2  It  j 

21.655  22.026  5280.3)55  0.04120  0.58  7.0 

15  15  7  26 

22.885  147.009  5280.8593  0.02927  0.53  7.0 

6  6  j  3  14 

24.008  34.427  ,5281.2275  0.02424  '0.J8  7.0 

8  7  1  8  |  j 

25.127  284.325  5281.5755  0.02223  1  0. 3-  }  7 . # 

3  5  2  7  1 

26.683  182.741  5282.0896  I  0.02080  I  0.50  7.0 

13  12  4  16  ! 

27.950  320.727  5282.3761  0.01746  }o.)6  7.7 

8  8  1  6  1 

29.175  101.303  5282. b085  0.03667  !0.58  5.9 

14  11  5  18 

30.235  3.401  5283.5021  0.05260  0.32  5.9 

9  4  2  9 

31.344  270.317  5284.1356  0.04248  |  0.63  6.9 

17  16  3  13  | 

32.197  296.754  5284.6377  0.04846  !  0.60  5.6 


33.181  326.465  3285.2688 

19  18  5 


0.05581  i  0.64  ]  6 . 3 

21  i 


725.0  Aug  5  7044.5810  0.040210  62.7379 
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10,  5516 
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31 

26 

j 

12 

7 
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7 46  .0 

7043.9807 

0.04041 1 

62.7402 

147.9)42 

39.285 

22  3.25  1 

3 

1  1 

8 

11 

16 

15 

754.0 

Sep  3 

7043. 7319 

0. 040486 

62.7414 

322.0416 

40.489 

42.709 

J  3 

2 

8 

6 

9 

12 

1  1 

761  .0 

10 

704 . .4281 

0.0405)2 

62.7380 

299.  )96)  1 

1  41.58  7 

336.842 

,  ^  1 

2 

9 

6 

1  J 

1  2 

|  76  7  .0 

16 

7043.0906 

0.040578 

62.73)8  1 

!  279,98  35 

1  m2. 52b 

25.47  ) 

I  l 

35  | 

2 

9 

6  j 

6  ■ 

1  2 

773.0 

22 

7042.6663 

0.040587 

62.73)4  | 

260. 5642  1 

4  3.428 

76.614 

:  36  ' 

1 

2 

8 

7  ! 

8 

1  3 

12 

7  79.0 

28 

7041 ,9341 

0.040579 

62. 7344 

241,1  )95 

44 , 386  | 

1  31  .473 

_ 1 

3 

10 

7  1 

7 

15  | 

14 

5285,9900 

4 

0.0431b  | 
29  i 

0.5  3 

1 

i  4.6 

!  »s 

5286.2633 

) 

0.02764  j 

18  j 

i  0 . h0 

!  4.9 

j  H4 

j 

5286,  b()74 

0.02497 

1  1 

0.49 

.  6  .M 

I  N\» 

528b  .909  3 

1  J 

0.01973  1 

9 

0.47 

■  H .  0 

N-, 

528  7.1)65  j 
3 

0.01  304 

16 

0.61 

6  .6 

Ml 

5287.28  )4 

3  | 

o .  :>  1 1  b  7 
ih 

0.6  3 

6.2 

'  M- 

5287.56  38  j 

0. 0 J  9h7 

9  \ 

0.-W 

■  7 . 9 

H> 

5287. ‘1056  ' 
l  l 

0.02^66 

1  5 

0.45 

5 , 9 

M  3 

5208.2855  , 

0.03115 

0.m7 

b  .  0 

:  ‘*0 

5288.76)4  i 

0.05089  ' 

0.47  ! 

5.0  , 

76 

J  semi  major  axis  (km) 
e  eccentricity 
i  inclination  (deg) 

'*  ri*ht  ascension  of  node  (deg) 

-  argument  of  perigee  (deg) 

*  orbits  with  Hewitt  camera  observations 


i  :  21  l 

85  u .08715  |  ' 

3  15  1 


mean  anomaly  at  epoch  (deg) 
mean  motion,  n  (deg/day) 

third  coefficient  in  polynomial  for  mean  anomaly 
measure  of  fit 

time  coverage  of  observations  (days) 
number  of  observations  used 
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2 . 3  Motion  of  perigee 

Since  the  inclination  is  close  to  63°,  the  perigee  moves  very  slowly,  at  less  than 
0.2  deg/day,  and  the  argument  of  perigee  -  increases  from  2°  to  44°  in  the  263  days 
between  first  and  last  epochs,  as  shown  in  Fig  l.  The  odd  zonal  harmonics  in  the  geo¬ 
potential  have  a  great  effect  on  eccentricity  when  i  =  62.74°:  the  amplitude  of  the 
oscillation  is  nearly  50  km,  so  a  decrease  in  perigee  distance  of  over  30  km  is  to  be 

expected  as  u  increases  from  2°  to  44°.  In  fact,  a(l  -  e)  decreases  from  6787  km  on 

the  first  orbit  to  6756  km  on  the  last.  This  corresponds  to  a  decrease  in  perigee 

height  over  a  spherical  Earth  from  409  km  initially  to  378  km  at  the  end.  Fig  1.  Over 
an  oblate  Earth,  the  corresponding  values  are  409  km  and  386  km. 

2 . 4  Observational  accuracy 

The  residuals  of  the  observations  on  the  first  32  orbits  are  summarized  in  Table  2 
for  stations  with  5  or  more  observations  accepted.  The  residual  of  an  observation  is 
a  combination  of  the  observational  error  and  any  error  in  the  orbit,  and  the  value  given 
in  the  Table,  the  rms ,  produces  a  bias  towards  the  larger  values,  which  for  visual  obser¬ 
vations  relative  to  the  stars  are  usually  observations  made  in  poor  conditions  of  seeing. 
For  these  reasons  the  capability  of  visual  observers  in  good  conditions  is  usually  reckoned 
to  be  about  half  the  rms  residual.  For  the  US  Navy  station  29,  the  angular  residuals  are 
geocentric,  and  need  to  be  multiplied  by  a  factor  of  about  5  for  comparison  with  the  other 
( topocentric)  observations.  All  observers  with  at  least  one  observation  accepted  have 
been  sent  copies  of  their  residuals. 

Table  2 

Residuals  for  observing  stations  with  more  than  5  observations  accepted 


Station 

Number  of 
observations 
accepted 

Rms  residuals 

Range 

km 

Minutes  of 

arc 

RA 

Dec 

Total 

1  US  Navy 

187 

1  .7 

1  .9 

2.6 

2  US  Navy 

130 

2.  1 

2.2 

3.0 

3  US  Navy 

1  34 

2.0 

1  .9 

2.8 

4  US  Navy 

141 

1 .6 

2.0 

2.6 

5  US  Navy 

161 

1  .9 

1  .  7 

2.6 

6  US  Navy 

175 

1  .6 

1.8 

2.4 

29  US  Navy 

625 

0.5 

0.  3* 

0.4* 

4  l  4X  Cape  Town 

28 

2.3 

2.  3 

3.2 

2122x  Malvern  5 

1  3 

1  .6 

1  .2 

2.1 

2  1  25x  Street 

7 

2.9 

1 

3.3 

2I55X  Bahrein  2 

8 

2.8 

4.4 

5.2 

2265x  Farnham 

6 

2.0 

l  .  3  1 

2.4 

2303  Malvern  Hewitt  camera 

9 

0.02  ! 

0.02 

0.0  3 

241 4X  Bournemouth 

59 

3.6 

3.5 

5.0 

2420x  Willowbrae 

28 

1  .** 

2.  1 

2.8 

2577  Cape  kine theodol i te 

75 

o.8  ; 

0.7 

1  .  ! 

4I68-1-  Vries 

1 1 

4.0 

3.8 

5.5 

6702x  Jokioinen 

2  1 

!  i.O 

|  3.6 

4.  ; 

*  Geocentric 
x  Visual  stations 


THE  29:2  RESONANCE 
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3 . I  Analysis  of  inclination,  i 

The  theory  for  29:2  resonance  is  detailed  in  Refs  3-5,  where  all  the  parameters 
used  here  are  defined.  The  first  term  of  the  inclination  equation  is 


di 

dt 


r— ~  )3°  ( 


(29  -  2  cos  i)F 


30,29,14  l  29 


0,2 


sin  $  +  C 


.0,2 

29 


•)  • 


a) 


where  $  =  2(ui  +  M)  +  29  (ft  -  v) 

is  the  resonance  angle,  v  being  the  sidereal  angle.  The  values  of  $  and  $  are  given 

in  Fig  2:  at  exact  resonance  b  *  0.  The  use  of  two  extra  terms  in  (i)  -  taking  (y,q)  ” 

(1,0),  (1,1)  and  (1,-1)  in  the  notation  of  Ref  4  -  gave  indeterminate  results,  very 

probably  because  the  ($  +  m)  terms  interfered  with  the  $  terms  as  a  result  of  the  very 

slow  variation  of  u>  .  The  fitting  of  i  was  therefore  made  with  equation  (1),  in 

2 

integrated  form,  using  the  THROE  computer  program  . 

Before  being  fitted  by  THROE,  the  36  values  of  inclination  from  Table  1  were  first 
cleared  of;  zonal  harmonic  and  luni solar  perturbations  (combined  maximum  value  being 
0.0070°),  using  the  PROD  program  with  numerical  integration  at  daily  intervals;  tesseral 
harmonic  perturbations,  determined  by  PROP  (maximum  0.0016°);  and  atmospheric-rotation 
perturbations  (maximum  0.0018°)  determined  within  THROE,  using  an  atmospheric  rotation 
rate  A  of  1.0,  which  gave  a  better  fitting  than  the  other  alternative  tried,  A  =  0.9. 
Earth  and  ocean  tide  effects  were  not  taken  into  account,  and,  in  recognition  of  this, 
the  sd  of  one  value  of  inclination,  on  orbit  4,  was  degraded  from  0.0003°  to  0.0005°. 

A  density  scale  height  H  of  60  km  was  used,  appropriate  to  a  mean  height  of  430  km 
(0.75  H  above  mean  perigee  height). 

The  36  modified  inclination  values  were  then  fitted  by  the  integrated  form  of 
equation  (!),  using  THROE.  In  the  first  fitting  the  measure  of  fit  e  was  1.5,  so 
the  nine  worst-fitting  values  were  degraded  successively,  two  by  a  factor  of  3  and 
seven  by  a  factor  of  2.  Fig  3  shows  the  final  fitting,  and  the  values  with  their 
original  standard  deviations,  the  nine  degraded  values  being  marked  by  horizontal  bars. 

The  values  of  the  lumped  coefficients  finally  obtained  were: 

0,2  02 

lo9C29  =  -  10  i  15  ,  1 °9S2g  -  -  76  ±  12  ,  (2) 

with  *  =  0.88.  These  were  within  1  sd  of  the  values  obtained  on  the  first  fitting. 

Kig  1  shows  that  the  overall  effect  of  the  resonance  was  to  increase  the  inclination  by 

about  0.003°. 

.0,2 

lor  1 9  7 l -b2E  the  numerical  expression  for  in  terms  of  the  individual 

coefficients,  as  obtained  from  the  RAE  computer  program  PROF,  is 


30,29 


15C 


32,29 


,4»CJ4i2,  .  °.82C)6i29  -  0.70C38i29 


0'6‘C40.29  * 


<3) 


c 

■sj 

o 


with  the  same  equation  for  S  ,  on  replacing  C  by  S  throughout. 
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)  .  2  Analysis  of  eccentricity,  e 

In  attempting  to  analyse  the  variations  in  eccentricity,  the  36  values  of  e  from 
Table  I  were  first  corrected  for  lunisolar  perturbations,  using  PROD,  and  then  the  values 

tMj)  j  -  CM 1 ) 

of  M.,  were  modified,  being  replaced  by  — — - —  - — j —  ,  where  (M  )  is  the  value 

*"  n  + 1  n  11 

of  on  the  nth  orbit,  at  epoch  .  This  allows  for  the  integrated  effect  of  drag 

between  successive  epochs^,  as  required  by  THROE.  (This  correction  is  not  significant 
in  fitting  i  .)  In  the  THROE  runs  the  scale  height  H  has  to  be  taken  at  a  height 
! . 5H  above  perigee,  and  a  value  of  65  km  was  used. 

The  appropriate  equation  for  fitting  the  variation  in  e  due  to  resonance  is 

/ 

( 


de 

dt 


n(— 

\a; 


29 


(JJ 

'  ,6F29.29,U)C29  Sin<' 


J  ,  J 

j)  -  S.^  cosd  -  u] 


+  I2F 


h"*’3  . 


29,29,1 3)^29 


sin($  +  cj)  -  S 


29 


cos  ( 


(4) 


taking  (y,q)  =  (1,1)  and  (1,-1).  Accordingly,  the  values  of  e  ,  after  removal  of  air- 
drag  and  zonal -harmonic  perturbations  within  THROE,  were  fitted  by  THROE  using  equation  (4) 
in  integrated  form.  Unfortunately  the  values  obtained  for  the  lumped  coefficients 
(C,S)^^  and  (C,S)”^’^  were  indeterminate ,  probably  as  a  result  of  being  highly 
correlated  because  ui  varies  so  slowly.  Also  the  value  of  c  was  high  (3.0),  and,  even 
after  five  values  of  e  had  been  degraded  in  accuracy,  thus  reducing  t  to  2.1,  the 
values  of  the  lumped  harmonics  were  still  all  less  than  twice  their  sd. 

The  first  possible  escape  route  from  the  'correlation  trap'  is  to  ask  whether  one 
pair  of  terms  in  (4)  is  small  and  can  be  ignored.  For  1976-62E,  the  PROF  computer  pro¬ 
gram  gives 
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J  ,  1 

"29 

_-l,3 
29 


C29,29  “  6‘  UC3I  ,29 


10#80Cn,29  '  4‘69C35,29  -  5*08C37,29 


*  C 


29,29 


* .  38 C 


3  f  ,29 


37C 


3  3,29 


1  *  1  8C  3  5 , 2  0  -  -4IC37,29 


(5) 


...  j 


and  similarly  for  S  .  Equations  (5)  indicate  that  the  most  probable  value  for  the 


-1,3/  1,1 


ratio  C 


29 


C  is  about  1,  and  similarly  for  S  .  Since  the  multiplying  factors 


outside  the  two  eurlv  brackets  in  equation  (4)  have  numerical  values  of  0.44  and  0.76 

. ; ,  I 

respectively,  the  most  probable  situation  is  that  the  terms  are  of  about  the  same 


magnitude  as  the  C 


terms  in  (4),  so  that  neither  can  legitimately  be  ignored. 


Another  escape  route  might  be  found  by  raking  the  magnitudes  of  the  two  terms  as 
equal.  Since  they  are  opposite  in  sign,  the  equation  (4)  would  reduce  to  the  form 
1,1  1,1 


K  <  (’ 


-  S  sin  ’)  sin  ,  where  K  is  constant.  So  if 


increased  from,  say. 


and  fit  a 


to  !I0°,  it  would  be  possible  to  take  a  constant  value  of  sin  ~ 
l  ,  ,q)  =  (1,0)  variation.  But  unfortunately  sin  u,  increases  steadily  from  0  to  0.7 
and  cannot  be  taken  as  constant.  Thus  w  is  too  nearly  constant  to  allow  separation 
vu  Lite  effects  of  the  lumped  coefficients,  but  not  constant  enough  to  allow  useful 
s  i  snp  l  i  f  i cat  ions . 

Even  if  one  of  these  escape  routes  had  been  open,  however,  there  would  still  be 
another  obstacle  to  face:  the  amplitude  of  the  oscillation  in  perigee  height  is  very  large 
for  1976-62E  (about  46  km),  and  even  the  latest  set  of  odd  zonal  harmonics^  may  well  have 
uncertainties  uf  up  to  2  km  at  this  inclination.  Since  the  change  due  to  resonance  is 
likely  to  be  less  than  I  km,  it  is  not  possible  to  remove  the  odd  zonal  harmonic  perturba¬ 
tion  with  adequate  accuracy. 

This  difficulty,  together  with  the  correlation  between  the  lumped  coefficients, 

prevents  a  satisfactory  solution.  But  the  values  of  eccentricity  still  need  to  be 

allotted  some  fitted  curve  for  future  comparison;  so  it  was  decided  to  adjust  the  values 

of  the  odd  zonal  harmonics  by  altering  J 7  so  that  (a)  the  value  of  e  was  near  minimal 

and  (b)  the  values  of  the  lumped  harmonics  were  consistent  with  the  maximum  credible 

-  -  “52. 

values,  found  by  inserting  (C,S)„  =  10  / £  in  (5)  and  taking  the  sum  of  the  numeri- 

cal  values  of  the  terms.  (The  individual  coefficients  are  generally  less  than  10  / 1“  - 
see,  for  example,  Ref  8  -  and  they  will  rarely  all  be  of  the  same  sign.)  This  gives 
max . ( C , S) ^  ~  300  x  10  ^  and  max. (C ,S) ^ ^  —  150  *  10  ^  .  The  best  solution,  which 

gave  values  within  1  sd  of  these  limits  and  had  a  reasonably  low  e  ,  namely  e  =  1.6, 

-9  -9 

was  obtained  by  changing  from  the  standard  value  in  THROE,  -326  x  10  ,  to  -306  *  10  . 

The  resulting  values  of  lumped  harmonics  were: 


9-1,1  9-1,1 

10*C29  =  340  ±  750  10  S29 

9_-l,3  q-H ,3 

10  C29  =  -  320  ±  430  10  S9 


(6) 


Fig  4  shows  the  values  of  e  ,  with  their  original  sd,  and  the  fitted  curve.  In  the 
fitting,  three  of  the  standard  deviations  were  relaxed  by  a  factor  of  3,  and  nine  by  a 
factor  of  2.  These  twelve  values  are  marked  by  horizontal  bars.  The  fitting  is  fairly 
sat i s far tory ,  and  it  is  probable  that  determinate  values  for  the  lumped  coefficients 
would  have  emerged  if  they  had  not  been  so  strongly  correlated  and  of  similar  magnitude. 

Since  the  variation  of  perigee  height  with  time  (Fig  1)  is  near-linear,  another  way 
uf  removing  any  error  in  the  calculated  amplitude  of  the  odd  zonal  harmonic  perturbation 
would  be  to  include  a  linear  term  in  the  fitting.  This  was  tried,  but  the  numerical 
•■ouf f icient  of  the  linear  term  was  less  than  its  sd,  and  the  general  appearance  of  the 
fitting  was  not  appreciably  altered. 

Finally,  the  first  five  values  of  e  ,  which  seem  rather  low,  were  omitted,  and 

i  fitting  with  only  the  second  pair  of  lumped  coefficients  was  made.  The  resulting 

-9 

'lines,  with  .1^  =  -  296  *  If)  ,  were: 
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9-"1’3  9--U3 

10  C29  =  -  462  *  43  10  S?9  =  302  '  46  ,  (7) 

with  r  =  1.39.  Tin*  standard  deviations  in  (7)  are  much  lower  than  in  (6),  but  cannot  be 

-9 

accepted  as  realistic  because  a  change  ot  J ^  to  -27b  x  ]Q  changes  the  two  values  to 
-528  and  -49  respectively. 

1  ^  The  most  useful  information  to  emerge  is  that,  in  all  the  21  fittings  attempted, 

C0Q  was  always  strongly  negative.  Since  its  numerical  value  seems  unlikely  to  exceed 

-9  -9 

150  y  10  ,  the  results  point  towards  a  value  somewhere  near  -150  *  10 

Although  the  present  analysis  of  e  has  not  yielded  good  values  of  lumped 
harmonics,  this  is  mainly  because  of  the  lack  of  knowledge  of  values  of  the  odd  zonal 
harmonics.  When  these  are  better  established,  a  successful  analysis  may  be  possible, 
since  the  orbital  data  are  basically  accurate  enough. 

4  COMPARISON  WITH  GEM  1 0B 

The  only  comprehensive  gravity  model  that  goes  to  order  and  degree  36  is  the 

9  .  .10.. 

Goddard  Earth  Model  1 0B  and  recent  tests  with  accurate  resonant  orbits  indicate 

-9 

that  the  15th-order  terms  in  GEM  1 0B  are  probably  accurate  to  3  x  !0  .  The  terms  of 

-9 

order  29  and  degree  30-36  are  likely  to  be  less  accurate,  and  an  error  of  5  *  10  may 
be  tentatively  assigned.  (The  use  of  GEM  10C,  which  is  the  same  as  GEM  10B  to  order  and 
degree  36  but  goes  to  degree  180,  is  not  very  useful  here,  because  the  terms  of 
degree  >36  do  not  greatly  influence  the  lumped  coefficients.) 

-9 

Assuming  an  error  of  5  *  10  and  ignoring  the  terms  of  degree  38,40,...,  GEM  1 0B 

gives 

9-0,2  9-0,2 

10  C99  =  -  8  ±  15  10  S29  =  -  1 1  -  1 5  (8) 

Comparison  with  the  values  (2)  from  1976-62E  shows  good  agreement  for  C  and  disagreement 

-0,2 

for  S  .  The  high  negative  value  of  S0g  obtained  from  1976-62E  derives  directly  from 

the  main  increase  in  inclination  between  MJD  43615  and  43675  (see  Fig  3),  where  270°  < 

v  <  320°  (see  Fig  2)  and  sin  0  in  equation  (1)  therefore  has  a  value  between 

-0.64  and  -1.0.  This  ’main  increase’  in  inclination  seems  securely  based,  because  it 

-0,2 

is  so  strongly  confirmed  by  all  the  outlying  points  in  Fig  3.  So  the  value  of  S^9 

from  equations  (2)  seems  preferable  to  that  from  equations  (8). 

-"1,3 

The  value  of  from  GEM  1 0B  is  not  likely  to  be  very  realistic,  because 

of  the  neglect  of  terms  of  order  37,39  ,...;  but  the  >.  =  33  term  in  GEM  10B  gives 

9-” 1 ,3 

a  contribution  of  -56  t  20  to  10  C^9  ,  which  is  consistent  with  the  strongly  negative 

value  indicated  by  1976-62E. 

5  CONCLUSIONS 

The  orbit  of  Cosmos  837  rocket  has  been  determined  at  36  epochs  spread  throughout 
the  first  nine  months  of  1978,  when  the  effects  of  29:2  resonance  were  being  felt.  About 
3100  observations  were  used,  of  which  about  757  were  supplied  by  the  US  Navy. 
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The  3  b  orbits  obtained  are  given  in  Table  1  and  show  chat  the  sd  in  inclination 
varies  between  0.0005°  and  0.0017°  while  the  sd  i"  eccentricity  varies  from  4  -  10  h  to 
I ^  x  10  .  The  best  standard  deviations  are  each  equivalent  to  30  m  in  position. 

The  36  values  of  inclination  (cleared  of  zonal-harmonic,  J0  9  ,  air-drag  and 
lunisoiar  perturbations)  were  fitted  with  a  least-squares  theoretical  curve  to  give  the 
following  values  of  lumped  29th-order  coefficients: 

9-0,2  9-0,2 

10  C.,9  -  -  10  t  15  10  S^9  =  -  76  •  12  . 

The  analysis  of  eccentricity  was  not  so  successful  because  of  (a)  the  difficulty 

of  accurately  removing  the  very  large  perturbation  (31  km)  due  to  odd  zonal  harmonics, 

and  (b)  the  interference  between  the  two  pairs  of  lumped  harmonics,  caused  by  the  very 

slow  variation  of  the  argument  of  perigee  at  this  inclination  (62.74°).  However,  there 

— — 1*3  —  9 

is  a  strong  indication  that  C^9  is  strongly  negative,  of  order  -150  *  10 

The  values  of  lumped  harmonics  from  1976-62E  -  which  appear  to  be  the  best  so  far 
obtained  at  29:2  resonance  -  will  be  used  in  future  determinations  of  individual 
coef f icients,  when  results  for  a  variety  of  inclinations  are  available. 
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These  will  be  used  in  future ,  when  enough  results  at  different  inclinations  have 
accumulated,  to  determine  individual  coefficients  of  order  29.  The  values  of 
lumped  harmonics  obtained  from  analysis  of  the  values  of  eccentricity  were  not 
well  defined,  because  of  the  high  correlations  between  them  and  the  errors  in 
removing  the  very  large  perturbation  (31  km)  due  to  odd  zonal  harmonics. 
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